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FOREWORD 


The  initial  investigations  of  plastic  models  described  in  this 
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Atomic  Support  Agency  (DASA).  The  current  work  on  variable  moduli  models 
is  part  of  Contract  No.  DACA  39-67-0-0048,  "Investigation  cf  Ground  Shock 
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Army  Engineer  Waterways  Experiment  Station  (WES)  under  DASA  sponsorship. 

This  work  has  been  done  in  close  cooperation  with  the  WES  Project 
Engineer,  Mr.  J.  G.  Jackson,  Jr.  The  authors  wish  to  thank  him  for  his 
assistance  in  describing  the  behavior  of  real  materials.  They  also  wish 
to  acknowledge  the  contribution  of  Dr.  Hans  H.  Bleich,  Consultant,  to 
the  early  theoretical  development  of  the  plastic  and  the  variable  moduli 
models . 

The  WES  contract  is  under  the  general  supervision  of  Mr.  W.  J. 
Turnbull,  Chief,  Soils  Division.  COL  Levi  A.  Brown,  CE,  is  the 
Contracting  Officer. 


ABSTRACT 


A  historical  outline  of  mathematical  models  used  previously 
to  represent  soil  and  rock  behavior  in  ground  shock  computations 
at  low  pressure  levels  is  presented.  The  problems  and  de¬ 
ficiencies  of  these  early  approaches  are  discussed.  Sub¬ 
sequently,  two  types  of  models  are  developed;  one  a  plastic 
model  in  which  the  yield  condition  depends  on  the  mean  stress 
and  in  which  different  variable  bulk  moduli  are  used  in  loading 
and  unloading.  The  second  type  of  model,  called  the  "variable 
moduli  model",  has  variable  shear  as  well  as  bulk  moduli,  but 
no  explicit  yield  condition.  The  behavior  of  both  types  of 
model  in  uniaxial  strain  and  triaxial  compression  tests  is 


examined . 
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INTRODUCTION. 


The  growing  interest  in  the  determination  of  the  response 
of  structures  at  relatively  high  pressure  contours  has  stimu¬ 
lated  the  study  of  ground  shock  from  both  direct  and  air  in¬ 
duced  sources.  As  the  investigations  have  progressed  into 
the  range  of  thousands  of  psi  air  overpressure  contours,  the 
studies  of  combined  air  and  direct  ground  shock,  as  well  as 
air-induced  ground  shock  only  have  required  the  development 
of  considerably  more  comnlicated  and,  hopefully,  more  accurate 
material  models.  'Jhereas  previously,  air-induced  ground  shock 
effects  at  low  overpressure  levels  were  approximated  by  means 
of  linearly  elastic  models,  the  requirements  of  problems  in¬ 
volving  combined  air  and  direct  ground  shock  effects  necessi¬ 
tated  a  definition  of  the  material  behavior  in  a  range  ex¬ 
tending  from  the  hydrodynamic  behavior  at  extremely  high 
pressure  regions  in  tne  neighborhood  of  Ground  Zero,  through 
various  types  of  inelastic  solid  material  regions  and  finally 
to  essentially  elastic  ranges  at  sufficiently  low  pressures. 
Parallel  to  the  theoretical  development  of  these  more  realistic 
material  models  have  been  the  difficult  demands  on  the  ex¬ 
perimentalists  for  tests  which  would  mirror  the  behavior  of 
the  material  throughout  the  entire  range  of  pressures.  The 
tests  served  two  major  purposes:  (a)  to  give  an  idea  of  the 

prooer  type  of  material  behavior  as  an  aid  in  the  mathematical 
modeling  of  the  soil  or  rock  at  appropriate  pressure  levels, 
and  (b)  to  provide  data  from  which  the  required  material 
constants  in  the  mathematical  models  could  be  determined. 


In  the  high  pressure  range  (say  above  100  KB)  Hugoniot  data 
is  obtained  from  high  velocity  impact  flying  plate  tests.  In 
the  low«r  range  from  about  5  KE  to  SO  KB,  data  is  presently 
being  obtained  from  both  flying  plate  tests  and  explosive 
tests  with  spherically  diffusing  wave  geometries.  Below  2  KB, 
the  usual  soil  tests  of  uniaxial  strain  and  triaxial  compression 
are  utilized.  Unsupported  uniaxial  compression  tests  which  are 
run  for  rock  are  essentially  a  special  case  of  a  triaxial  com¬ 
pression  test  in  which  the  lateral  stress  0^  is  zero. 

For  a  given  soil,  the  uniaxial  strain  and  triaxial  com¬ 
pression  tests  give  the  most  detailed  information  and  are 
generally  utilized  in  the  development  of  the  mathematical  model, 
particularly  in  those  cases  in  which  air-induced  ground  shock 
only  is  being  studied.  Typical  experimental  curves  for  each 
test  are  shown  in  Fig.  (1).  The  stress-strain  curve  in  the 
uniaxial  strain  case,  Fig.  (la),  typically  shows  a  reversal  in 
curvature  on  loading.  On  unloading,  the  slope  is  almost 
constant  and  is  much  larger  than  the  slope  during  initial 
leading,  except  for  a  sharp  tail  in  the  low  stress  range. 
Reloading  generally  follows  the  unloading  curve  up  to  the 
previous  maximum  stress  and  then  continues  along  the  initial 
loading  curve.  The  lateral  stress  0^  >  requited  to  maintain 
uniaxial  strain,  is  sometimes  also  measured,  Fig.  (lb). 

In  the  triaxial  compression  tests,  the  stress-strain 
curves,  Fig.  (Ic),  are  concave  downward  up  to  a  horizontal 


__  3  -- 


tangent  at  what  is  called  "failure".  Different  curves  are 
obtained  for  different  values  of  the  confining  stress  . 
Although  the  lateral  strain  is  not  generally  available, 

experimental  work,  currently  in  progress,  is  aimed  at  providing 
this  information.  By  plotting  the  Mohr's  circles  at  failure 
for  different  lateral  stresses  cr^  ,  the  Mohr  envelope,  Fig.  (Id), 
is  obtained.  This  envelope  is  generally  either  a  straight  line 
(for  some  dry  sands)  or  is  concave  downward  (for  partially 
saturated  soils)  .  Both  the  theoretical  and  the  experimental 
programs  have  made  considerable  progress  in  recent  years. 

The  present  paper  will  describe  some  of  the  recent  theoretical 
approaches  which  have  been  and  are  being  used  in  the  development, 
of  the  material  models. 

Before  discussing  these  models,  it  is  perhaps  of  some 
interest  to  review  the  development  of  material  models  from  a 
historical  viewpoint.  Originally,  linearly  elastic  models 
were  used  in  the  analysis  of  air-induced  ground  shock  effects 
at  low  pressure  levels  (20-100  psi  contours).  As  technological 
changes  in  weapons  design  and  delivery  systems  forced  the 
locations  of  the  hardened  structures  to  considerably  higher 
overpressure  levels,  inelastic  effects  were  introduced  into 
the  analysis.  For  the  most  part,  these  effects  consisted  of 
introducing  a  yield  condition  into  the  material  model  of  the 
form  suggested  by  Prager  and  Drucker,  Ref.  [ 1  ] : 


aJl  +  V  3 2  =  k 


(1-1) 
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or  a  condition  of  the  von  Mises  type 

*  k  (1-2) 

where  a  is  a  coefficient  which  depends  upon  the  angle  of  in¬ 
ternal  friction  of  the  material,  k  is  a  constant  pronor tioned 
to  the  cohesion,  =  -3p  is  the  first  invariant  of  the 
stress  tensor,  and  is  the  second  invariant  of  the  stress 
deviators.  The  early  models  of  materials  which  were  assumed 
to  act  linearly  elastic  below  certain  stress  configurations 
and  to  obey  the  yield  condition  of  Eq.  (1-1)  were  called 
Coulomb  materials  with  cohesion  or  Prager-Drucker  materials. 

A  considerable  amount  of  both  numerical  and  theoretical 
work  has  been  done  for  the  special  case  in  which  the  cohesion 
constant  k  was  taken  to  be  zero,  Refs.  [2],  [3],  (41,  and  (5). 

Such  materials  have  been  called  Coulomb  materials  and  satisfy 
the  yield  condition 


For  cases  in  which  the  yield  condition  satisfies  Eq.(I-l)  or 
Eq.(l-3),  the  materials  have  a  conical  yield  surface  and  ex¬ 
hibit  the  phenomenon  of  dilation,  that  is,  an  inelastic 
volume  increase  when  loaded  to  sufficiently  high  stress  in 
compression  and  subsequently  unloaded.  Figure  (2)  shows 
typical  stress-strain  plots  for  such  materials  in  a  uniaxial 
strain  compression  test  configuration.  It  is  seen  that  upon 
full  unloading,  a  net  residual  extension  is  obtained  for  this 


t 
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essentially  one-dimensional  configuration.  In  three-dimensional 
problems  (two  dimensional  space-time),  this  manifests  itself 
in  an  inelastic  increase  of  the  volume.  The  von  Mises  material, 
on  the  other  hand,  which  satisfies  the  yield  condition  of 
Eq .  (1-2),  has  a  cylindrical  yield  surface  and  exhibits  no 
dilatancy  effects  as  shown  in  Fig.  (2)-  Theoretical  solutions 
to  several  problems  have  also  been  obtained  for  the  von  Mises 
material,  Refs.  [6],  [7],  and  [8J.  Problems  arose  in  the  use 
of  both  of  these  models.  First  and  perhaps  foremost  was  the 


dilatancy  effect  from  the  dependence  of 


/  J  2  '  on  Jj^ 


in  the 


yield  condition  for  the  Coulomb  material  model.  While  such 
phenomena  have  been  observed  in  certain  types  of  sands  under 
special  loading  conditions,  this  is  certainly  not  typical  of 
many  of  the  soils  and  rocks  which  are  encountered  in  practical 
design  problems.  The  von  Mises  material  on  the  other  hand, 
requires  that  the  material  maintain  essentially  the  same  shear 


strength  at  all  stress  levels,  since  no  variation  in  \j  J ^  is 
allowed  in  the  yield  condition  as  the  stress  range  increases 
and  hence,  jj^j  increases.  While  such  a  situation  is  probably 
quite  true  at  extremely  high  pressure  levels,  e.g.,  as  the 
material  passes  from  a  hydrodynamic  to  a  solid  state,  it  is 
certainly  not  generally  true  at  lower  stress  levels  in  which 


tests  usually  show  that  an  increase 


ll 

in  /  J .  before  yieli 


occurs  as  increases.  Consequently,  two  problems  arose: 
(1)  while  many  materials  of  interest  acted  essentially  as  a 


Coulomb  material  at  low  pressures,  at  higher  pressures  they 
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acted  as  a  von  Mises  material  and  consequently,  the  mathe¬ 
matical  material  model  would  have  to  act  in  a  similar  fashion, 
and  (2)  the  Coulomb  model  would  have  to  be  amended  so  that 

it  had  a  zero  net  dilation  at  the  end  of  a  complete  loading- 
*) 

unloading  cycle 

The  next  advance  in  model  development  was  the  intro¬ 
duction  of  a  nonlinear  pressure-volume  relation,  i.e.,  a 
variable  bulk  modulus  K  =  KCJ^),  in  the  material  description. 
At  the  same  time,  the  shear  modulus  G  was  kept  constant  ' . 

The  use  of  the  variable  bulk  modulus,  in  addition  to  giving 
the  proper  concavity  to  the  stress-strain  curve  from  a  uni¬ 
axial  con pression  test,  aided  materially  in  obtaining  a  smooth 
transition  at  high  pressures  from  the  hydrodynamic  material 
(governed  by  a  suitable  equation  of  state)  to  the  solid  in¬ 
elastic  material.  In  the  neighborhood  of  the  transition,  the 
solid  inelastic  material  acts  as  a  von  Mises  materiil  with  a 
nonlinear  pressure-volume  relation.  By  properly  choosing  the 


In  any  material  in  which  the  plastic  flow  rule  is  derivable 
from  a  yield  condition  which  depends  on  the  mean  pressure  p 
(or  the  first  invariant  J^),  the  incremental  dilation  can¬ 
not  be  zero  but  the  net  dilation  from  a  complete  loaditg- 
unloading  cycle  may  be  eliminated. 

A  model  in  which  the  shear  modulus  G  was  made  a  function 
of  in  unloading  was  also  considered,  but  extensive 
numerical  ground  shock  calculations  were  not  made  for  this 


model . 
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constants  in  Che  expression  K  ■  K(J^),  the  shock  velocities 
for  the  inelastic  solid  material  in  the  selected  transition 
range  may  be  matched  to  those  given  by  the  equation  of  state 
of  the  hydrodynamic  material  in  this  range,  thus  eliminating 
possible  impedance  mismatches  and  resulting  in  a  imooth  tran~ 
sJ  tion. 


The  problem  of  modeling  a  material  which  obeyed  a  Coulomb 
type  yield  condition  at  low  pressures  and  a  von  Misci  yield 
condition  at  higher  pressures  was  solved  by  requiring  a  smooth 
transition  between  the  two  conditions  through  the  use  of  a 
yield  condition  of  the  form 


f(Jfj)  +  y  J,  -  k 


U-4) 


where  the  function  of  was  chosen  so  that  the  material  had 
the  proper  behavior  in  the  appropriate  ranges,  e.g,,  es¬ 
sentially  Coulomb  at  low  values  of  the  mean  pre,sure  p  and 
von  Mises  at  high  values.  In  an  early  version  of  this  type 
of  model  for  a  volcanic  tuff,  the  relation 
J, 


«(1  +  ~~)2  Jj  ^ 


»  k 


r — •-» 

1  J2  B  k 


Jlf  <  C 


J1  >  c 


(1-5) 

(1-6) 


was  utilised,  Fig.  (3),  together  with  the  relations 


K  -  Kq  -  Vl  +  (1-7) 

C  ®  G 

o 


(1-8) 
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In  this  model  the  constants  K  ,  G  ,  a.  and  k  were  chosen 

o  o 

from  whatever  information  was  available  from  triaxial  com¬ 
pression  (a  and  k)  and  uniaxial  strain  (K  and  G  )  tests- 

o  o 

The  constants  K^  and  Kj  were  chosen  to  make  the  shock  veloci¬ 
ties  from  the  hydrodynamic  equation  of  state  model  (which  was 
used  in  the  higher  pressure  ranges)  match  those  from  the 
solid  inelastic  material  model  over  a  20  KB  to  100  KB  range 
of  mean  pressures.  From  Fig.  (3)  it  is  seen  that  the  yield 
condition  in  this  model  was  chosen  such  that  a  smooth  tran¬ 
sition  exists  between  the  Coulomb  type  material  and  the 
von  Mises  type  material.  Since  the  cohesion  constant  k  was 
kept  the  same  at  both  low  and  high  values  of  J^, a  reversal 
in  curvature,  as  shown,  was  obtained  for  the  yield  surface. 

This  model  was  inherently  questionable  because  of  this  nega¬ 
tive  curvature  of  the  yield  surface,  since  uniqueness  has 
not  been  mathematically  proved  for  such  cases.  It  represented 
one  of  th*2  very  early  taodels  of  this  type  and  has  long  been 
superseded  by  what  we  consider  to  be  more  accurate  models 
y’lich  are  described  in  detail  later  in  this  report. 

The  later  models  utilized  a  yield  condition  of  the  form 

-/T7  ->i 

|/  J2  «  k  -  “did  +  2^)  Jj  +  c  >  0  (1-9) 

«  k  +  --■  J1  +  c<0  (1-10) 

as  shown  in  Fig.  (4).  Again,  the  condition  approximates  the 
Coulomb  material  at  low  pressures  and  a  v»  a  Mises  material 
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at  high  pressures,  while  the  transition  surface  is  such  that 
the  curvature  is  positive  over  the  entire  surface,  a  sufficient 
condition  for  uniqueness.  A  detailed  discussion  of  a  model  in 
which  the  yield  condition  is  utilized  is  given  in  Section  II 
for  an  alluvium-playa  material. 

The  use  of  different  nonlinear  pressure-volume  relations, 
i.e.,  nonlinear  bulk  moduli  in  loading  and  unloading  enables 
one  to  eliminate  the  net  dilation  for  a  complete  loading¬ 
unloading  cycle  in  uniaxial  strain,  Fig.  (5).  This  also 
satisfies  the  requirement  of  a  steep  unloading  slope  as  found 
in  the  experiments. 

It  is  in  order  to  make  some  general  observations  regarding 
the  problems  involving  dilation,  particularly  those  which  arise 
in  reloading  situations.  In  material  models  which  utilize 
plastic  yield  conditions,  it  is  usually  advisable  to  derive 
the  flow  rule  from  a  plastic  potential  which  is  identical  to 
the  yield  condition,  since  for  such  cases,  uniqueness  theorems 
have  been  derived  for  specific  yield  conditions  such  as  the 
von  Mises  and  Coulomb  conditions,  and  can  be  deduced  for  other 
more  generalized  cases.  The  introduction  into  the  model  of  a 
nonlinear  pressure  -  volume  relation  together  with  the  yield 
condition  creates  no  new  difficulties  when  the  same  curve  is 
used  for  both  loading  and  unloading.  While  this  may  be  per¬ 
missible  in  some  ranges  for  certain  types  of  rocklike  materials, 
there  are  many  other  rock  and  soil  materials  for  which,  in  the 
range  of  in'srest,  material  property  tests  indicate  important 
hysteretic  energy  losses  in  a  complete  loading-unloading  cycle. 
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In  mathematical  models  for  such  materials,  different  K  and  G 

m 

relations  must  be  used  in  loading  and  unloading.  This  may 
lead,  however,  to  problems  on  reloading  as  shown  in  Fig.  (6). 

For  example,  if  we  load,  unload,  and  then  reload  along  the  un¬ 
loading  curve  until  the  original  loading  curve  is  reached,  we 
will  generally  find  ourselves  off  the  original  loading  curve 
because  of  the  inherent  dilatancy  effect  (inelastic  volume 
change)  in  the  material  if  a  yield  condition  is  used  in  which 
the  first  invariant  of  the  stress  tensor  appears.  The 
question  of  how  to  return  to  the  original  loading  curve  arises 
and  is  not  an  easily  answered  one.  For  example,  one  way  of 
doing  this  would  be  to  use  a  yield  condition  of  the  form 
f  (J  )  +  "I/  J  =>  k,  but  at  the  same  time  use  a  flow  rule  derivable 

1  1  2  .r~~r 

from  a  plastic  potential  of  the  type  F  *»  u  -  k.  Such  a 

procedure  would  still  determine  yielding  by  the  general  yield 

condition,  but  would  correct  for  the  yielding  by  changing  only 

the  stress  deviators,  and  no  change  would  be  made  in  the  mean 
J1 

pressure  p  =  -  —  in  the  correction  for  yielding.  This  procedure 
would  eliminate  any  incremental  dilation  effect  in  any  portion 
(time  step)  of  the  numerical  calculations.  This  approach  has 
been  used  in  several  of  the  existing  codes.  Another  possibility 
might  be  to  return  to  the  original  F-V  curve  by  some  sort  of 
interpolation  procedure.  A  third  possibility  would  be  to  con¬ 
tinue  loading  along  a  curve  parallel  to  the  original  P-V  curve. 

There  are  several  possibilities  no  one  of  which  can  be  said  to 
be  correct  as  compared  to  the  physical  evidence  for  real  materials. 
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It  is  important,  however,  to  determine  the  sensitivity  of  the 
various  models  to  the  several  possible  procedures.  Investigations 
of  these  effects  are  currently  being  made. 

An  alternative  approach  to  the  determination  of  mathematical 
models  has  led  to  a  series  of  what  we  call  "variable  moduli" 
models.  In  these  cases,  the  bulk  modulus  K  and  the  shear 
modulus  G  are  assumed  to  be  functions  of  appropriate  stress 
and/or  strain  tensor  invariants, and  no  plastic  yield  conditions 
are  specified  for  the  material.  Different  functional  relation¬ 
ships  for  K  and  G  are  used  in  loading  and  unloading.  The 
study  of  such  models  is  relatively  recent  and  has  been  con¬ 
fined  at  this  time  to  theoretical  investigations.  It  is  hoped 
that  the  work  will  be  applied  to  the  analysis  of  real  materials 
in  ground  shock  problems  in  the  near  future. 

The  models  are  developed  by  requiring  that  the  truncal ed 
series  expansions  for  both  K  and  G  in  terms  of  stress  and/or 
strain  tensor  invariants  have  their  constants  chosen  so  that 
the  three-dimensional  model  conforms  to  both  the  triaxial  com¬ 
pression  and  the  uniaxial  strain  results  when  it  is  run  for 
these  configurations.  Once  the  appropriate  form  of  the  K  and 
G  functions  is  determined,  the  constants  appearing  in  these 
expressions  are  determined  by  means  of  a  curve-fitting  routine 
from  uniaxial  strain  and  triaxial  compression  test  data  for  a 
given  soil. 
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One  possible  advantage  of  these  models  is  that  they 
accurately  reproduce  the  available  test  data  from  the  material 
property  tasts  that  are  generally  made.  This  is  an  advantage 
over  the  previously  discussed  plastic  models  in  which  K  was  a 
function  of  and  G  was  kept  constant.  In  the  latter  case, 
excellent  reproductions  of  uniaxial  strain  results  were  ob¬ 
tained,  but  the  results  from  triaxial  compression  tests  were 
not  reproduced.  The  combination  of  the  two  types  of  models, 
namely  a  variable  modulus  model  which  also  satisfies  a  yield 
condition, has  not  yet  been  studied  in  detail. 

Section  III  describes  three  of  the  variable  moduli  models 
that  have  been  studied  thus  far.  Conclusions  are  presented 
regarding  their  suitability  as  models  for  real  soil  materials. 
While  the  studies  on  the  variable  moduli  models  are  still  at  a 
relatively  early  stage,  it  appears  that  this  approach  shows 
promise . 
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II  NONLINEAR  PLASTIC  MATERIAL  MODEL  -  ALLUVIUM- PLAY A • 


The  model  utilizes  a  nonlinear  bulk  modulus  and  a 
constant  shear  modulus,  together  with  a  modified  Prager- 
Drucker  yield  condition.  The  vield  condition  is  of  the  form 


k  -  aJi  (i  +  =  k 


a(Ji)Ji 


J1  +  c  >  0  (II-l) 


=  k  + 


"effective 


Jx  +  c  <  0  ( 1 1-2) 


This  condition  approximates  the  Pragor-Dr ucker  condition 

for  <<  c,  while  in  the  high  pressure  range,  the 

material  behaves  as  a  von  Mises  material  with  k  =  k  c c  , 

effective 

Fig,  (4).  At  =  -c,  a  smooth  transition  occurs.  It  is 

of  interest  to  note  that  in  this  model,  the  value  of  the 

cohesion  constant  k  at  low  pressures  is  less  than  the  corre- 

etc 

spending  constant  ^effective  =  k  +  ~  at  higher  pressures, 
thus  resulting  in  a  transition  yield  surface  which  maintains 
a  positive  curvature,  a  sufficient  condition  for  uniqueness. 


From  available  test  information,  In  particular  uniaxial 
strain  test  data,  it  was  noted  that  this  material  exhibited 
a  considerable  amount  of  hysteresis  in  a  loading-unloading 
cycle.  In  order  to  obtain  the  proper  amount  of  hysteresis 
in  the  mathematical  model,  different  pressure-volume  relations 
were  used  for  loading  and  unloading.  These  relations  were 
defined  by  variable  bulk  moduli  through  the  relation 


J,  ■  JKWj)  ik„ 


(11-3) 


The  shear  modulus  G  =  G  was  kept  constant  for  initial 

o 

loading . 


In  uniaxial  strain,  experimental  information  showed 
that  this  material  is  characterized  bv  an  initial  softening 
and  a  subsequent  hardening,  Fig.  (la).  The  simple  ouadratic 
relation  for  the  bulk  modulus,  Eq .  (1-7),  used  previously 
in  the  earlier  tuff  model,  was  not  sufficient  to  both  match 
the  hydrodynamic  material  shock  velocities  at  high  pressures 
and  to  provide  the  proper  curvature  for  the  stress-strain 
curve  at  lower  pressures.  Consequently,  an  altered  bulk 
modulus  was  used  for  higher  pressures  | J ^ |  >  Yc*  At  lower 
pressures,  | J ^  J  <  yc,  the  simplest  expression  which  yields 
an  initial  negative  curvature  in  ur'axial  strain  and  which 
matches  the  value  of  the  higher  pressure  K  relation  at 
|jjJ  =  yc  was  used.  The  initial  loading  bulk  modulus  K  is 
given  by  the  relation 


■WV  ■  [Ko  -  f  Kl>  +  y!  IJ1  -  ¥>: 


<  yc  (11-4) 


=  Kq  -  Kl(Jx  +  yc)  +  k2(Jx  +  yc)2  >  yc  (II-5) 


so  that  K,  _  *»  K  =  elastic  bulk  modulus  when  J,  **  0  and  again 
LD  o  I 

when  ■  -yc .  Between  these  values  of  drops  to  a 

minimum  at  J1  =  -  ^  .  At  *  -yc,  both  and  its 

^  ^LD 

derivative  with  respect  to  Jj  ,  i.e.,  —  ,  are  continuous. 


Above  =  -yc,  Kld  increases  monotonica 1 ly . 
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At  very  large  pressures,  one  might  expect  the  material 
to  behave  essentially  as  a  fluid  (little  or  no  energy  dis¬ 
sipation)  with  similar  loading  and  unloading  characteristics. 
A  quantity  J.  Tn  is  defined  as  a  measure  of  the  pressure 
above  which  the  behavior  of  the  material  in  loading  and  un¬ 
loading  is  similar.  Below  J.  n  ,  particularly  where  ex- 
perimental  data  is  available,  it  appears  that  the  slope  of 
the  unloading  curve  in  uniaxial  strain  is  considerably 
greater  than  the  corresponding  loading  curve  slope,  thus 
indicating  a  considerable  amount  of  energy  dissipation  in 
this  material.  Moreover,  the  unloading  slope  appears  to 
diminish  gradually.  The  bulk  modulus  for  unloading  was 
chosen  to  give  this  behavior,  and  is  given  by 


WJ1>  ■  Ko  -  K1J1  -  Vi 

J1  +  J1,LIQ  <  0 

( 1 1  —  6  ) 

’  K4  Y  1  -  t  J1 

J1  +  J1 , LIQ  >  0 

( I X  —  7  ) 

where 


K  +  J.  ...  K.  +  j:  ...  K„ 
o _ l.LIQ  1 _ 1,LIQ  2 


LIQ 


( 1 1 — 8  ) 


is  chosen  so  that  K„,,  is  continuous  at  J,  =  -J,  ,  Tr.  . 

UN  Jv  1  1 ,  LIQ 

41  iJN 

iiOcc  uiai.  tilt  (juaitbii.)'  ^  j  x  o  u  ioCOuLiiiUOus  au  u  u  x  »■»  fiuxiib  » 

k 

The  value  of  K„„  can  never  be  zero  since  J,  <  —  bv  the 

u  w  i  a 

yield  condition  requirements.  When  |Jl|  >  J l.LIQ  >>  YC* 

Ky,{  is  always  greater  than  K^p  ,  but  the  difference  is  small. 
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Below  =  - J ^  very  gradually  decreases  with 

decreasing  pressure  until  at  very  low  pressure  de¬ 

creases  rapidly.  Figure  (7)  shows  a  typical  diagram  of  the 
relation  between  K  and  for  both  loading  and  unloading. 

With  respect  to  the  shear  modulus  G  in  unloading,  three 
different  variations  of  the  model  have  been  considered. 

First,  the  same  constant  Gq  was  used  in  both  loading  and  un¬ 
loading  for  a  series  of  ground  shock  calculations  for  an 
alluvium-playa  medium.  Secondly,  a  constant  value  of 
G  =  GyN  which  differs  from  the  initial  loading  shear  modulus  Gq 
was  studied  theoretically.  Finally,  a  variation  in  which 
the  shear  modulus  on  unloading  has  been  allowed  to  vary  as 
a  function  of  over  a  portion  of  the  region, 

G  =  G  when  s .  ,  A  e .  .  >  0  ( 1 1  —  9 ) 

o  ij  i  j 

G  =  (1  +  — —  ]  when  S;lj  A  eij  <  0  (II-1Q) 

has  been  studied  theoretically. 

The  determination  of  the  various  constants  which  appear 
in  the  model  was  based  predominantly  on  uniaxial  strain  test 
data.  Attempts  to  make  the  model  also  fit  available  in¬ 
formation  on  the  Mohr  failure  envelope  from  triaxial  com¬ 
pression  tests  were  unsuccessful.  Since  the  geometry  for 
the  ground  shock  problems  are  closer  to  a  uniaxial  strain 
test  configuration  than  to  that  of  a  triaxial  compression 
test,  the  model  was  chosen  to  conform  as  closely  as  possible 
to  the  uniaxial  strain  experimental  data.  The  constants  K 
and  G^  were  chosen  from  the  uniaxial  strain  data,  while  the 
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constants  K^  and  were  determined  so  that  in  tlv  higher 
pressure  ranges  (say  20  KB  -  90  KB)  the  shock  velocities 
of  the  solid  material  matched  those  from  the  h- drodynamic 
equation  of  state  material  over  this  range  of  pressures. 

The  constant  y  was  chosen  so  that  the  strain  at  the  point 
of  inflection  of  the  uniaxial  stress-strai.i  cur^e  falls 
within  the  range  of  the  experimental  valves.  The  quantity  c 
was  determined  from  available  strength  data  while  the 
constant  J.  .  T_  was  essentially  guesseJ.  As  new  experimental 
results  into  higher  pressure  ranges  become  available  under 
the  present  DASA  experimental  program,  it  is  hoped  that 
the  J,  T T_  value  may  be  determined  directly  from  the  test 
data.  The  plastic  constants  a  and  k  were  approximately 
determined  from  triaxial  compression  data  and  were  modified 
to  make  the  model  conform  to  the  uniaxial  strain  results. 


Theoretical  curves  for  th<  uniaxial  strain  test  are 

presented  in  Figs.  ( 8 )  —  ( 9 )  for  the  plastic  model  described 

in  this  section.  The  results  are  for  the  case  in  which  G„„ 

UN 

is  given  by  Eqs.  C 1 1  •“  9 )  -  (II  10).  The  following  constants 
were  used  in  the  calculations 


G  =  0.5  KB 
o 

K  =  1.222  KB 
o 

K1  =  8.0 
K2  =  0.07 (KB) 


Onlv  two  of  the 
independent . 


0  =  (  . 15 
k  =  0.001  KB 
c  =  2  5  KB 

*) 

(VQ=  0.32) 


four  constants,  K  , 
o 


3 

Po  =  1.6  gm/cm 
y  =0.15 

J 1 , LIQ  ”  12  KB 
(E0=  1.3200  KB) 

G  .  v  ,  and  E  are 
o*o  o 
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Figure  (8)  shows  the  curve  of  the  axial  strt.es  a^,  versus 
the  axial  strain  for  a  complete  loading-unloa ling  cycle 
for  the  uniaxial  strain  configuration.  Figure  (9)  shows 
the  lateral  stress  versus  the  axial  stress  0^  in  uniaxial 
strain.  It  is  seen  that  both  curves  g  iva  the  characteristic 
shape  of  the  experimental  results,  Figs,  (la)-(lb). 

A  series  of  ground  shock  calculations  using  the  three- 

dimensional  (two-dimensional  space-time)  alluviurr-p]  ava  model 

have  been  made  for  a  model  configuration  in  which  (a)  the 

same  constant  shear  modulus  G  =  G  has  been  used  in  both 

o 

loading  and  unloading,  and  (b)  the  yield  condition  of 
Kq.  (II-l)  has  been  used,  but  a  plastic  potential  of  the 
form  fA  -  k  =  0  has  been  used  In  computing  the  plastic 
strains,  so  as  to  avoid  the  dilation  problems  which  occur 
on  loading  and  unloading  as  described  in  some  detail  in 
Section  I. 

The  effect  of  using  G  =  G^  in  both  loading  and  unloading 
results  in  a  stress-strain  curve  for  the  uniaxial  test  which 
although  slightly  softer  in  unloading,  is  essentially  the 
same  as  the  curve  in  Fig.  (8).  The  lateral  stress  curve  for 
unloading  would  differ  from  the  curve  shown  in  Fig.  (9),  since 
the  unloading  and  loading  curves  would  coincide  at  low  stresses. 
The  use  of  the  generalized  yield  condition,  Eq.  ( 1 1 - 1 ) ,  together 
with  a  plastic  potential  which  is  independent  of  ,  results 
in  a  slightly  softer  stress-strain  curve  than  is  shown  in 
Fig.  (3). 
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It  is  of  interest  to  consider  the  behavior  of  plastic 
material  models  of  this  general  type  in  the  triaxial  com¬ 
pression  cest.  Figure  (10)  shows  a  sketch  of  a  typical 
ci^  -  0^  versus  diagram.  Because  of  the  variable  bulk 
modulus  K(p),  the  initial  loading  case  first  softens  and 
then  hardens.  As  the  yield  surface  is  reached  (at  different 
levels  for  different  values  of  the  lateral  stress  0^)>  the 
material  "fails"  and  "flows"  along  the  horizontal  lines 
shown  parallel  to  the  strain  axis.  It  is  apparent  that 
plastic  models  of  this  type  in  which  a  constant  shear  modulus 
G  is  used  do  not  mirror  the  usual  stress-strain  curves  from  a 
triaxial  compression  test. 

A  way  to  remedy  this  situation  is  an  approach  in  which 
variable  moduli  are  used  for  bath  K  and  G.  This  approach 
is  discussed  in  detail  for  several  cases  ir  Section  HI. 
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III  VARIABLE  MODULI  MATERIAL  MODELS. 


This  section  describes  the  development  of  the  variable 
moduli  models  in  which  no  explicit  plastic  yield  condition  is 
specified.  As  mentioned  in  the  Introduction,  these  models 
are  an  alternate  to  the  approach  utilizing  plastic  yield  con¬ 
ditions  as  described  in  Sections  I  and  II.  In  the  variable 
moduli  models,  both  the  shear  modulus  and  the  bulk  modulus 
will  be  assumed  to  depend  upon  the  stress  and/or  strain  in¬ 
variants  . 


The  mathematical  description  of  the  model  is  in  terms 
of  the  incremental  stress-strain  relations 


s1}  -  2G  k±j 


p  =  3K  e 


(III-l) 


(III-2) 


where  and  e^  are  the  deviatoric  stress  and  strain, 

* ) 

respectively,  and  p  and  e  are  the  mean  stress  and  strain. 
In  general,  different  functions  G  and  K  apply  in  initial 
loading  and  subsequent  unloading  and  reloading.  The  present 
discussion  will  be  largely  confined  to  the  case  of  initial 
load ing . 


It  is  noted  tha:,  even  in  cases  of  initial  loading, 
there  is  not  in  general  a  unique  stress-strain  relation. 

The  final  state  of  strain  depends  not  only  upon  the  final 
state  of  stress,  but  also  upon  the  stress  path  used  to  reach 


In  this  section  compressive  stress  and  strain  are  defined 
as  positive  in  accordance  with  the  usual  soil  mechanics 
convention . 
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the  final  state.  In  this  sense,  the  variable  moduli  material 
cannot  be  considered  a  nonlinear  elastic  material  where  such 
a  unique  stress-strain  relation  would  exist.  Three  different 
models  or  this  type  are  discussed  in  this  section:  (1)  a 
model  in  which  Poisson's  ratio  is  kept  constant;  (2)  a 
model  in  which  both  K  and  G  are  given  as  a  function  of  the 
strain  tensor  invariants,  and  (3)  an  essentially  equivalent 
model  to  (2)  in  which  K  is  a  function  of  the  strain  tensor 
invariants  and  G  is  a  function  of  the  stress  tensor  invariants. 
The  latter  model  is  more  convenient  for  large  computer  codes 
since  it  requires  less  storage  capacity  per  computation  point 
than  does  the  model  of  type  (2). 

Each  of  these  models  will  be  discussed  from  the  stand¬ 
point  of  the  uniaxial  strain  and  triaxial  compression  tests, 
since  these  are  the  soil  tests  which  will  generally  be  available 
for  the  determination  of  the  material  constants  for  the  three- 
dimensional  (two-dimensional  space-time)  codes.  It  should 
be  noted  that  when  such  models  are  run  for  the  simple  geome¬ 
tries  of  these  two  material  tests,  it  is  a  necessary  condition 
that  the  experimental  test  results  be  closely  checked  by  the 
mathematical  computations.  It  will  be  shown  that  the  constant 
Poisson's  ratio  model  of  type  (1)  cannot  match  both  the  uni¬ 
axial  strain  and  the  triaxial  compression  experimental  results, 
but  that  both  of  the  models  of  type  (2)  or  type  (3)  do  match 
both  tests,  and  consequently,  appear  to  be  suitable  mathematical 


material  models. 
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( 1 )  Constant  Poisson's  Ratio  Model. 

A  simple  variable  moduli  model  which  has  been  suggested 
and  used  by  some  investigators  is  one  in  which  the  ratio  ^ 
remains  constant.  It  will  be  shown  here  that  a  model  of 
this  type  cannot  satisfy  both  the  uniaxial  strain  and  the 
triaxial  compression  test  results. 

By  analogy  with  the  linearly  elastic  relation 

K  2 ( 1+v) 

G  ~  3 (l~2v) 

this  model  may  be  called  a  constant  Poisson's  ratio  material 

However,  K  and  G  will  not  be  constants,  but  functions 
of  the  mean  strain  e 


K  «  K(e) 

G  -  G  (e )  =  |  K  (e  ) 

where  Ms  a  constant.  Alternatively,  K  and  G  may  be 
viewed  as  functions  of  the  pressure  p.  Since  for  initial 
loading,  the  pressure 

e 

p  =  |  3K(0  d£ 

0 

is  a  monotonic  function  of  e,  a  unique  inverse  e(p)  exists. 
Thus,  the  moduli  K  and  G  may  also  be  written  as 

K  -  K(p)  =  K[e(p>] 

G  =  G(p)  =  £  K{ e (p) ] 


(III-3) 


(III-4) 


<  1 1 1  -  5 ) 


( 1 1 1  —  6 ) 
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(1 1 1-14) 
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Integrating  Eq.  (III-12)  and  Eq.  (Ill— 14) ,  the  expressions 
for  the  principal  stresses  and  a ^  become: 

el 

°l(el)  =  (1  +  |  Kt3)  dK  (III- 15) 

0 

£1 

03UL)  =  (1  -  |)  j  K(|)  d£  (III-16) 

0 

Equations  (III-12),  (III-15),  and  (III-16)  uefine  a 
constant  Poisson's  ratio  material  during  initial  loading 
in  a  uniaxial  strain  configuration.  If  K(e)  initially 
decreases  and  subsequently  increases,  the  axial  stress  - 
axial  strain  curve  [see  Eq.  (Ill— 12) ]  will  have  the 
characteristic  sof teni.ig-hardening  shape  found  experi¬ 
mentally,  Fig.  (la). 

From  Eqs.  (III-15)  and  (III-16), 

°3(el)  =  (|?f)  °l(ei)  ( I I I- 17 ) 

so  that  the  radial  stress  versus  axial  stress  relation  is  a 
straight  line  of  slope  •  This  agrees  reasonably  well 

with  experimental  data,  Fig.  (lb). 


(1 11-13) 


Thus,  =  2p  and  =  3p.  The  axial  strain  increment  is 
given  by 


O  -  • 

G  =  — i.  +  £_ 

1  2G  3K 


1  11 
( —  +  — 1  — 
VG  3K;  3 


(IIt-19) 


and  the  slope  of  the  stress-strain  curve  by, 

do. 


de,  "  6+1 


K(e) 


(111-20) 


The  radial  strain  increment  is 


e  a  -  —  +  £—  =  ( _  _i_  _L\  _! 

3  4G  3K  1  2G  3K;  3 


(111—21) 


or 


£  .  -.(g.-2,)  •  =  _v£ 

e3  2(6+1)  L 1  V61 


(111-22) 


which,  not  surprisingly,  is  the  same  as  the  linear  elastic 
relation.  The  mean  strain  increment  is  also  proportional 
to  e . 


1  » 


e  =  3  *£1  +  2£3 ' 


-L.  £ 

6+1  1 


(111-23) 


so  that  the  actual  total  mean  strata 


:n  by 


e  =  eo  +  6+1  Ael 


(m-24) 


where  e^  is  th.e  strain  reached  during  the  loading  of  the 

specimen  hydrostatically  up  to  the  stress  o^iSnd  where 

Ae,  =  e,  -  e  is  the  measured  axial  strain. 

1  1  o 


Typical  experimental  triaxial  results,  Fig.  lie), 


have  the  following  characteristics: 
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1.  The  stress-strain  curve  at  a  given  value  of  0^  is 
concave  downward. 

2.  At  some  point  a  horizontal  tangent,  "failure" ,  is 
reached . 

3.  At  a  higher  value  of  a ^  the  initial  slope  and  the 

stress  difference  at  failure, (a,  -  a.)  ,  both 

1  3  max  ’ 

increase . 


From  Eqs.  (III-20)  and  (III-24)  one  sees  that  the  re- 

2 

d  O  i  jj. 

quirement  of  negative  curvature  - 2~  <  0  means  that  — 

de  ^ 

must  be  negative,  which  contradicts  the  behavior  of  K(e) 

which  was  required  for  the  uniaxial  test.  Moreover,  the 

assumption  K  >  0  is  not  compatible  with  the  existence  of 
da. 

"failure"  (-—  »  0). 


Even  if  a  yield  condition  of  some  type  were  included 

in  the  model  specif icaticn ,  well  below  the  yield  point, 

d  K 

negative  curvature  would  require  <  0.  On  the  other  hand, 
requiring  the  initial  slope  to  increase  with  ,  or 


de  jda-j 


9 

dK 

de 

3 

dK 

=  6-H 

de 

dp 

"  (B+1)K 

de 

Ae.=0 

J 

e 

p-d. 

- 

(HI-25) 


requires  that  ^  >  0,  wnich  contradicts  the  requirement 
that  ~  <  0. 


One  may  therefore  conclude  that,  although  the  constant 
Poisson's  ratio  model  appears  promising  when  only  the  uniaxial 
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results  are  considered,  such  a  model  must  be  rejected,  even 
in  conjunction  with  a  yield  condition,  since  its  behavior 
in  the  triaxial  configuration  contradicts  experimental 
evidence . 


( 2 )  Moduli  Dependent  Upon  Strain  Invariants . 

Models  in  which  the  moduli  were  assumed  to  be  independent 
functions  of  the  first  and  second  invariants  of  the  strain 
were  investigated  in  the  order  of  increasing  complexity. 
Initially,  strain  rather  than  stress  was  chosen  for  analytic 
convenience,  at  least  with  the  simpler  models. 


The  most  promising  of  the  strain  models  thus  far  con¬ 
sidered  is  defined  by  the  relations 


K  =  K(e) 


K  +  K.e  +  K„e 
o  1  Z 


G  =  G(e,  /i^) 


G  +  G. 
o  1 


I2  +  G  j e 


(III-26) 
(III  —  27) 


where  is  the  second  invariant  of  the  strain  deviators. 
Equations  (III-26)  and  (III-27)  may  be  thought  of  as  the 
first  terms  in  the  series  expansions  of  more  general 
analytic  functions  K  and  G  of  the  strain  invariants.  The 
quantity  |f  I2  was  used,  rather  than  I2  itself,  since  it 
is  of  the  same  order  as  e  and  the  components  of  the  strain 
tensor.  At  zero  strain,  the  bulk  and  shear  moduli  reduce 
respectively  to  Kq  and  Gq  ,  the  "linear  elastic"  values, 
which  are  related  in  terms  of  the  "elastic"  Poisson's  ratio  v 
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K  2 ( 1+v  ) 
o  _ o 

G  "  3 ( l-2v  ) 
o  o 

although  the  ratio  ^  Is  not,  in  general , constant .  With 
positive  and  G^  negative,  the  material  hardens  in  shear 
with  increasing  hydrostatic  strain  and  softens  with  in¬ 
creasing  shear  strain. 


(111-28) 


(a)  Uniaxial  Strain. 


The  general  incremental  stress-strain  relation  is 
given  by 


C*  S3 

eij 


JJL 


p6 


il 


2  [  G  +  G.  /iT  +  G.e]  3  ( K  +  K.e  +  K.e2] 
O  I  [  2  1  Ol  2 


In  the  case  of  uniaxial  strain,  =  S3  =  -  —  ,  so  that 


-s . 


e3  = 


4[Go  +  GL  yTJ1  +  G^e  ]  3[Kq  +  +  K2e2] 

fh  * 


■2  2  "1 


(since  in  ltad^g  e^  >  0)  and, 


el  =  2e  =  3  G1 


Combining  Eqs.  (III-30)  -  (III  —  32)  one  obtains  an  ex¬ 
pression  for  the  tangent  modulus  tt  any  strain 
K,  K, 

=  [  K  + 

1 


da 

de 


(111-29) 


0  (I 11-30) 


(III  —  31) 


( 1 1 1-32) 


K  K 

-  -  [Eo  +  -J  Ex  +  e2]  +  |  [Go  +  y  (Gx  +  /3  Gj) 


-  K  +  -  G  =  M 


(1 11-33) 
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(III  —  39) 
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and  since  j  I2  -  -y  (during  loading  e^  >  0),  the  strain 
rate  is  given  by 

3[K  +  K  e  +  K,e2) 

- - = -  £  (II 1-40) 

1Go  +  Gie  +  f  Vl> 

The  total  axial  strain  increment  is  given  by 


+  e 


(III  —  41) 


Combini*-  ,  Eqs.  (111-39)  and  (III-40),  the  tangent  modulus 
is  obtained  as 


or 


9[Kq  +  '  r  K2e2)[GQ  +  C,l&  +  'Ll 
3(Ko  +  Kxe  +  K0e2]  +  [Gq  +  G^  +  ■—  G  L  e  1 ) 


(111-42) 


do,  9K(e)  G(e,e.) 

dej  3K(e)  +  G  ( e  ,  e  L )  "  “ve,el' 


(III-43) 


where  E  may  be  considered  the  local  modulus  of  elasticity, 
i.e.,  Young's  modulus. 


Equations  (III— 40)  -  (III  —  42)  consitute  a  system  of 
nonlinear  d 1 r t i ren t ia 1  equations  relating  the  variables 
e,  ei:  and  .  The  system  may  be  solved  numerically 

by  forward  integration.  The  initial  conditions  are  obtained 
from  the  hydrostatic  state  of  stress  which  exists  in  a 
triaxial  compression  test  before  the  additional  axial 
loading  is  applied,  namely, 


J1 


°1(0)  =  °3 

e} (0)  =  0  (II 1-44) 

e , (0)  =  e(0)  =  e 
l  o 

where  ,  one  third  of  the  hydrostatic  volumetric  strain, 
is  the  sma 1  *  positive  root  of  the  cubic  equation 

0,  =  3K  e  +|  K .  e 2  +  K..e3  (111-45) 

3  oo  2  1  o  2  o 

Since  K  is  always  positive,  increases  monotonica lly 
wi  th  e 

o 

( c )  Choice  of  Constants  . 


If  all  stress  quantities  are  ncndimensionalized  with 
respect  to  the  initial  bulk  modulus,  ,  then  five  param¬ 
eters  remain  to  fully  describe  the  model,  namely, 


(III  —  46) 


The  first  of  these  is  inherently  positive  and  is  related  to 

the  initial  Poisson's  ratio  v  bv  Eq .  (Ill  —  28)  .  The  higher 

°  kl  K2 

order  terms  in  the  bulk  modulus  and  77—  may  be  positive 

K  K. 

0  O 

or  negative.  However,  the  values  are  restricted  by  the 


conu i t ion 
ii”'  .e 


minimum  K  occurs  at  e  . 

min 


Eq.  (111—26)  one 
,  which  is  posi¬ 


tive  (i.e.,  relevant)  only  for  and  K,,  of  opposite  sign. 


If  i’2  is  taken  to  be  negative, 


of  e  at  which  K  would  reach  a  maximum. 


would  be  the  value 
With  negative 
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and  positive,  t li e  requirement  that  the  minimum  bulk 
modulus  be  positive 


K 


K 


1 


min 


K  (e  )  =  K  -  7~—  >  0 
min  o  4K.2 


(III  —  47) 


requires  that 


K2  51  4K 


(III  —  48) 


K! 

It  should  be  noted  that  any  values  of  —  and  r are 

K 

0  O 

permitted  as  long  as  K  >  0  for  the  range  of  e  of  interest, 


especially  0  <  e  <  j 


In  the  triaxial  stress  test  the  initial  slope 


de, 


9 [ K  +  K.e  +  K0e  ) [G  +  G.e  ] 
c  lo  2  o  o  lo 

- - 2 - 

.  3 [ K  +  K.e  +  K.e  ]  +  (G  +  G.e  ] 

e^O  o  1  o  2  o  o  1  o 

e  =  e 


9K(e  )  G(e  ,0) 
_ o _ o _ 

3K (e  )  +  G(e  ,0) 
o  o 


(111-49) 


should  increase  with  increasing  lateral  stress  (or  e^) 
This  leads  to  the  relation 


3 [ K (e  ) ] 2  G.  +  [ G ( e  , 0 ) ] 2 ( K .  +  2K,e  )  >  0 

Ol  O  1  £  O 


(II 1-50) 


which,  for  e  sufficiently  small,  becomes 
o 


G1  1  Go.2  K1 

IT  +  3  (r>  r  >  0 

O  0  0 


With  <  0,  Eq .  (1 11  —  51)  sets  a  lower  bound  for 


(III  —  51) 
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K1  G1 

value  —  <  0  is  chosen,  —  is  taken  to  satisfy  Eq .  (Ill  —  51), 
o  o 

G1  K2 

and  jr-  to  satisfy  Eq .  (III-53).  Finally,  —  is  chosen  to 

o  o 

satisfy  both  Eqs.  (Ill  —  48)  and  (III— 56)  .  It  is  seen  that  the 

salient  features  of  both  the  uniaxial  strain  compression  test 

and  the  triaxial  compression  test  can  be  matched  by  material 

models  of  this  type. 


(d)  Numerical  Results. 


Typical  numerical  results  were  obtained  for  the  param- 
^1  ^2  G1 

eters  V  =  0.30,  —  =  -100,  -r—  =*  4000,  —  =  4  8  and 

O  K.  K  K 

O  O  O 

Gi 

=  -32.  For  the  uniaxial  corcpr  ssion  test,  the  plot  of 
o 

axial  stress  versus  axial  strain,  Fig.  (11),  shows  the  re¬ 
versal  in  curvature  which  is  characteristic  of  the  experi¬ 
mental  curves.  The  point  of  inflection  occurs  at  4.1% 

strain,  again  a  typical  experimental  value  The  plot  of  the 
radial  stress  required  to  maintain  uniaxial  strain  versus 
axial  stress  is  shown  in  Fig.  (12).  The  general  trend  of 
this  curve  conforms  to  the  experimental  curve,  Fig.  (lb). 

It  may  be  noted  that  a  slight  reverse  curvature  occurs  at 
the  low  end  of  the  theoretical  curve.  Although  a  change  in 
the  material  parameters  could  alter  or  possibly  remove  this 
effect  at  the  low  end,  it  is  not  completely  clear  as  to 
whether  such  effects  occur  in  some  real  materials.  In  any 
case,  this  effect  is  a  minor  one  and  is  not  sufficient 


grounds  for  rejecting  the  model.  Figure  (13)  is  a  plot  of 

the  deviatc-r  ~  versus  the  pressure  t?-  ,  again  for  uniaxial 
K  K. 

O  O 


strain.  Once  more,  at  the  low  end  there  is  a  reverse  in 
curvature.  However,  beyond  the  inflection  point,  the  slope 
of  the  curve  is  remarkably  similar  to  the  curve  of  the  playa- 
tyne  yield  condition,  described  in  Fig.  (4).  More  will  be 
said  about  the  similarities  and  differences  between  the 
plastic  and  variable  moduli  models  later. 

In  the  criaxial  test,  the  stress  difference  -  a ^  versus 
the  axial  strain  is  shown  in  Fig.  (14)  for  the  same  con¬ 
stants  that  were  used  in  the  uniaxial  calculation,  and  for 

°3 

—  =  0.04,  0.06,  0.08,  respectively.  The  curves  are  found 
o 

by  the  numerical  integration  of  Eqs.  (1II-4Q)  -  (III-42)  and 
reflect  the  experimental  results,  Fig.  (lc).  The  curves  are 
concave  downward,  exhibit  horizontal  tangents,  "failure",  and 
at  higher  lateral  stresses  the  initial  slope  increases  and 
the  stress  difference  at  "failure"  increases. 

Since  an  analytical  expression  for  the  stress  difference 
at  failure  is  unattainable,  Che  computer  program  was  run  to 
failure  for  nine  different  values  of  the  lateral  stress,  . 
Table  I  shows  the  stress  difference,  the  pressure, and  the 
axial  strain  at  failure  as  well  as  the  initial  modulus  for 
each  value  of  lateral  stress.  The  Mohr  envelooe  is  dotted 
in  Fig.  (15).  Again,  there  is  overall  agreement  with  ex¬ 
perimental  envelopes,  Fig.  (Id),  except  for  a  small  dis¬ 
crepancy  at  very  low  stresses. 


TABLE  I 


Triaxial  Test  at  Various  Lateral  Stresses 


(V 


o 


0.30, 


-100, 


K 

o 


4000, 


K 

o 


-32) 


No. 

°3 

(al“a3)failure 

^failure 

£  l 

f ai lure 

^initial 

K 

0 

K 

0 

K 

0 

k 

0 

1 

0.000 

0.0139 

0.0046 

0.0211 

1.200 

2 

0.001 

0.0153 

0.0061 

0.0232 

1.231 

3 

0.005 

0.0230 

0.0126 

0.0317 

1.350 

4 

0.010 

0.0577 

0.0292 

0.0667 

1.479 

5 

0.020 

0.1608 

0.0736 

0.1045 

1.594 

S 

0.040 

0.2385 

0.1195 

0.1212 

2.977 

7 

0.060 

0.2819 

0.1540 

0.1310 

3.973 

8 

0.080 

0.316? 

0.1854 

0.1374 

4.625 

9 

0.100 

0.3456 

0.2152 

0.1426 

5.117 
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On  the  basis  of  the  present  investigation,  it  appears 
that  the  variable  moduli  model  under  consideration,  defined 
by  Eqs.  (III-26)  and  (III-27),  shows  promise  as  a  mathematical 
reoresentation  of  the  mechanical  behavior  of  soils.  The  model 
behavior  in  both  uniaxial  strain  and  triaxial  stress  con¬ 
figurations  (in  initial  loading)  is  essentially  similar  to 
that  found  experimentally. 

No  attempt  was  made  to  use  this  model  to  match  actual 
data  numerically,  nor  was  the  problem  of  unloading  given  more 
than  a  cursory  glance.  These  activities  were  postponed  in 
order  to  develop  a  variable  moduli  model  which  depends  upon 

J  *) 

the  stress  invariants  p  and  j|  J2  as  well  as  the  mean  strain  e 
This  was  done  for  practical  as  well  a  theoretical  reasons. 
Practically,  in  the  general  three-dimensional  problem,  one 
is  physically  more  interested  in  stresses  than  strains  and 
already  stores  the  stresses  in  the  computer.  To  require  the 
storage  of  even  a  single  extra  quantity  at  each  grid  point 
significantly  reduces  the  number  of  grid  points  which  can  be 
utilized . 

Prom  a  theoretical  viewpoirt,  this  stress  invariant 
model  allows  an  interesting  comparison  with  plastic  models 
which  will  be  discussed  Idler. 


*) 


In  the  various  computer  codes  for  evaluating  ground  shock, 
the  density  p  is  generally  stored  at  each  point  and  the 
mean  stiair.  e  may  be  computed  from  it. 


(  3 )  Comb i n  ed  Stress-Strain  Variable  M  o  d  e 1 1  Mode  1 . 


This  section  presents  the  development  of  a  material 
model  in  which  the  bulk  modulus  K  is  a  function  of  the 
mean  strain  e,  and  the  shear  modulus  G  is  a  function  of 
che  «tress  tensor  invariants,  through  the  mean  pressure  n 
and  the  sauare  root  of  the  second  invariant  of  the  stress 


deviator 


fA- 


The  simplest  relation  in  which  the  shear  modulus 
deper.ds  upon  both  p  and  j 'j  is  given  by 

,I”- 

The  bulk  modulus  is  retained  as  a  quadratic  function  of 
the  volumetric  strain 

K  =  K  +  K.e  +  K„e2  (fee  Eq .  III-26) 

O  l  Z 

As  mentioned  in  Section  1 1 1  —  1  for  the  constant  Poisson's 
ratio  model,  a  unique  pressure-volumetric  strain  relation 
exists  during  initial  loading  as  long  as  the  bulk  modulus 
depends  only  upon  o  or  e.  Therefore,  Eq .  (111-26)  is 

equivalent  to  writing  K  as  a  function  of  p  [See  Eq.  (Ill  —  6)  J  . 

The  bulk  modulus  was  chosen  to  be  a  quadratic  in  c 
rather  than  in  d,  since  the  typical  uniaxial  strain  test 
curve  suggests  that  the  axial  stress  is  a  cubic  function 
of  the  axial  strain.  Since  the  moduli  K  and  G  refer  to 


the  incremental  stress-strain  relations, 
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p  =  3K  e  (See  Eq  .  1 1 1  - 1 ) 

®i1  =  2G  ®ij  (See  Eq‘  III_2> 

the  pressure  may  be  obtained  by  direct  integration  of 
Eq.  (1II-1)  as 

e 

D  =  3K(C)  d5  »  3KQe  +  |  K  ^  e 2  +  K2e3  (III-5B) 

0 

Moreover,  in  the  preparation  of  a  computer  code  for  the 
solution  of  three-dimensional  (two-dimensional  space-time) 
ground  shock  problems,  the  density  (or  the  volumetric 
strain  e^)  raust  be  stored  at  each  grid  point,  since  the 
density  appears  in  the  equations  of  motion.  Since  the 
relation  befaen  density  and  volumetric  strain  is  given  by 

E.  ,  «  3e  =  Hn  (-£-)  (III-59) 

kk  p 

Ci 

it  is  noted  that  the  use  or  a  bulk  modulus  which  is  a 
function  of  the  mean  strain,  Eq .  (Ill  —  26)  ,  requires  no 

A) 

significant  additional  computer  storage  .  It  should  also 
be  noted  that  in  Eq .  (Ill  —  57) ,  the  nondimensional  constant  y^ 
will  be  positive  for  a  G  which  hardens  with  increasing 
pressure,  while  y^  will  be  negative  for  a  G  which  softens 
an  Incroase  In  the  rfpviatorif  c tresses. 


In  the  large  ground  shock  problems,  the  problem  of 
minimizing  computer  storage  requirements  at  each 
grid  point  becomes  extremely  important,  since  the 
size  of  the  problems  challenges  the  capacity  of 
even  the  largest  computers. 


40  -- 


( a )  Uniaxial  S  t  r a i n , 


In  a  manner  quite  analogous  to  the  case  of  the  strain 
invariant  model.  Section  HI-2,  it  nav  he  shown  that  in 
uniaxial  strain 


=  3K  +  40  =  3  [  K'o  f  K  ^e  +  K2e2l  +  4(1',^  +  +  y  1 J  ,  )  (III-60) 

v  1  /T  /T 

Noting  that  ■  J ,,  =  s^  =  ~  (0j  -  p),  substitution  of 
Eq.  (111-58)  into  Eq.  (Ill- 60)  yields  the  first  order  non- 
homogeneous  differential  equation 

do  i  /a" 

~  2/1  Y1  Ol  -  (3Kq  +  4fio)  -  3[4Ko(Yi  -  y.)  +  Kj )  e  + 


+  3  ( 2K  x  (y  J  -  ~  y  x )  +  :<2  ]  e2  +  4K,{y,  -  ■—  y,)  (III-61) 


,  /I  -  ,  3 

'2 '  1  2  Y1)e 


Using  the  initial  condition  that  stress  and  strain  vanish 
simultaneously,  the  solution  for  stress  as  an  explicit 
function  of  strain  is  found  hv  integrating  Eq .  (III-61) 


20  y. 

+  -4  (K  + 

/3Yi  Y2  ° 


2/3  Yj  (2/3  Yx) 


- - y]>  [1  -  exp (2/3  v^el]  - 


2  /T  K 


i  f) 


/3  -  Y1  K2 

( y  -  V  )  +  — i  ( V  + - - - ) 

W1  2  Yl;  -2  '  1  rr  r. 


i  >  i  Ki  ,  /?  -  Yi 

|  ^1  “  7  V  +  -2  f' 

i  Yi  Y' 


n  1  ^Yli 

J  I  a  2f7  /g  - 

x  K  /  !*■  - - —  (Y,  -  ~r  Y,  )  e 

2  2J  /3Yj  1  2  1 


( I  11-62) 
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Closed  form  expressions  for  the  remaining  stress  quantities, 
Sj  and  0^  ,  are  easily  obtainable  from  Eqs.  (III-58)  and 


(111-62). 


S1  = 


E . 


2K 


2G  Y,  ... 

- P„  +  -1  [K  + - - -  + - z - 

[/3Yj  °  2/3  Y  (2/3  y  ) 


2 ]  >  [l  -  exp(2/3  Y ^e ) ] 


2/3  y,  Y,  K 

_  [ - 1  K  +  (r  +  -_/—))  e  _ 

-  O  -  L  1  —  - 


1 


1 


/3  Yj 


-  [-—  K  +  — j  K  ]  e 2  -  K  e3 

n  1  2  2 


(111-63) 


C  Yi  Ki 

d,  -  \  — — - —  +  — ~  (  K  +  - — - -  + 


2!<, 


3  1  /3  Yx  2yJ  °  2/T  Y 1  (2/3  Y ! ) 


-2-]  Ml  -  exp  (2/3  Y )  ]  + 


/3  K  Y,  K 

+  [—■—  (y !  +  /J  Yt>  +  (Ej  +  -r~-)l  6  + 


Y  - 

X 

/J  K, 


2Y 


/3  Y: 


Y  K 

+  [ - -1-  (Y,  +  /3  Y.)  +  -4  K  )  e2  +  —2—  (Y  +  /3  Y.) 

2Y 1  1  1  2y\  2  /)  Vj  1 


(1 1 1-64) 


Finally,  from  Eo  .  (HI-62),  the  slope  of  the  stress-strain 

curve  at  any  point,  or  the  tanpent  (constrained)  modulus 


is  found  to  be 
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--  4  4  -- 

stress  a  ^  .  It  should  also  be  noted  that  the  local  slope 
of  the  triaxial  stress-strain  curve 

^  0 1  9 kg  *) 

-  tkTg  2  E  (See  E<-  ul-43) 

goes  jo  zero  elien  0  -►  0,  so  that  the  stress  difference 

(a,  -  0,)  represents  a  point  of  horizontal  tangeney, 

1  3  max  r  ° 

i  .  e  .  ,  "failure". 

The  measured  strain  &e^  is  simnly  related  to  the 
strain  deviator  e^  ,  the  mean  strain  e,  and  the  initial 
(hydrostatic)  mean  strain  e^  by 

Ae,  =  e.  +  e  -  e  (111-73) 

Li  O 

where  e  is  found  for  the  given  lateral  stress  o,  by 
o  3 

Eq.  (Ill  —  43)  .  Equations  (111-57),  (III-26),  (  1  1  1-69), 

(III  —  43),  (III  —  73),  and  the  small  positive  root  e  of  the 
cubic  Eq .  (III-58)  completely  define  the  system  in  tri¬ 

axial  stress  for  all  valid  stress  states  o ^  ,  0 ^ 

l(0i  -  V  <  <0i  -  V™.,*1- 

If  the  Mohr  failure  envelope  were  plotted  fo  this 
material,  it  is  evident  from  Eq.  (Ill— 72)  that  the  plot 
would  be  a  straight  line  passing  above  the  origin,  similar 
to  the  yield  condition  for  a  Prager-Drucker  material.  This 
similarity  will  be  discussed  later  in  the  paper. 


It  can  be  shown  that  Eq.  (111-43)  applies  for  quite 
general  functions  K  and  O. 
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( c )  Choice  of  Constants . 

As  for  the  case  of  the  strain  invariant  model.  Section 
III-2,  five  parameters  are  required  to  fully  describe  the 
material  in  initial  loading: 


K. 


K, 


(III  —  74) 


The  choice  of  a  suitable  value  of  v  determines  the  ratio 

o 

G 

rr-  ,  Eq .  (111-23).  The  requirement  that  K  >  0  again  leads 
o 

to  the  condition  between  and  as  given  by  Eq .  (111-48). 
The  fact  that  y^  >  0  and  y^  <  0  for  physical  reasons  has 
already  been  discjssed,  as  has  the  inequality  between  them, 
Eq.  (III-70). 

To  further  restrict  the  range  of  the  five  material 
parameters,  one  requires  the  initial  slope  in  triaxial  com¬ 
pression  to  increase  with  lateral  stress,  and  the  initial 
curvature  in  uniaxial  strain  to  be  negative.  Differenti¬ 
ating  Eq .  (III-43)  with  respect  to  o ^ 


1  li_ 

9  do. 


3k2  35- 

do  . 


+  G 


2  sUL 

do  . 


( I  II -  7  5 ) 


(  U'  *  >  ' ) ' 


and  requiring  the  result  to  be  positive  when  0^  =  0^  , 
yields  the  inequality 


l 
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1Ko  +  Vo  +  ^2eo  ^  + 


+  [G^  +  Y,(3K  e  +  ~  K . e 2  +  K„e3)]2(K.  +  2K„e  )  >  0  (III-76) 

O  iOOidiO  i  0  I  Z  o 


For  Eq.  (Ill  —  76)  to  hold  in  the  limit  as  approaches  zero 
the  inequality 


Y1 


> 


( 1 1 1- 77  ) 


must  be  satisfied.  Equation  (III-77)  is  a  necessary  con¬ 
dition  for  the  initial  slope  in  the  triaxial  test  to  in¬ 
crease  with  increasing  lateral  stress  0^  .  Requiring  the 
initial  curvature  in  the  uniaxial  strain  test  to  be  negative 
results  in  the  condition 


K 

o 


< 


/3 

*1 


(III-78) 


Equation  (III  —  73)  is  obtained  by  evaluating  the  derivative 
of  Eq .  (III-65)  at  e  =  0.  an  attempt  to  find  an  analytic 
expression  for  the  inflection  point  leads  to  a  transcen¬ 
dental  equation  and  will  not  be  discussed  further. 


( d )  Unload  ing . 

The  model  for  unloading  (and  subsequent  reloading)  of 
variable  moduli  materials  is  presently  in  the  very  early 
stage  of  development.  Thus  far,  only  unloading  in  the  uni¬ 
axial  strain  configuration  has  been  considered.  Experi¬ 
mentally,  Fig.  la,  the  uniaxial  strain  unloading  curve  has 
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a  slope  much  larger  than  the  loading  slope,  which  is  ap¬ 
proximately  constant  until  very  low  stress  levels  are 
reached  . 

As  a  first  approach,  the  material  model  in  unloading 
was  chosen  to  be  defined  by  the  following  expressions  for 
the  bulk  and  shear  moduli: 

K  =  =  constant  (III-79) 

UN 

‘•‘OH’',*'!'*'.  K"(V  (UI-80) 

•  * 

where  h(J2)  is  the  unit  step  function.  The  effect  of  the 
step  function  in  Eq.  (111-80)  is  that  at  the  same  values 

I 

of  p  and  J2  ,  the  material  is  stiffer  in  shear  when  it  is 

•  t 

unloading  in  shear,  J2  <  0,  than  when  it  is  loading  in 

•  t 

shear,  J2  >0.  In  a  completely  general  three-dimensional 
configuration  the  terms  "loading"  and  "unloading"  no  longer 

have  such  clear-cut  meanings.  It  is  possible  that  the 

•  * 

material  will  be  loading  in  shear  (J2  >  0)  and  unloading  in 
pressure  (p  <  0)  simultaneously.  In  fact,  if  one  studies 
the  curve  for  as  simple  a  geometry  as  a  uniaxial 

strain  test,  Fig.  (lb),  it  can  be  shown  that  on  unloading, 
the  deviator  s^  =  y  (a^  -  0^)  [originally  positive]  first 
decreases,  then  changes  sign  apd  continues  to  decrease 
[increases  negatively]  until  a  minimum  is  reached  at  a  very 
low  stress  level.  Beyond  this  point,  appears  to  increase 
slightly,  i.e.,  to  decrease  in  absolute  value.  At  the  same 


—  48 


a,  +  2a3 

time,  che  pressure  p  =  - - -  and  the  strain  both 

decrease  monotonically  during  unloading.  Clearly,  the 
sharp  tails  found  experimentally  upon  unloading  at  very 
lot,*  stress  levels  in  both  the  uniaxial  stress-strain  curve, 

Fig.  (la),  and  the  radial  stress-axial  stress  curve,  Fig.  (lb), 
are  related  to  this  unusual  behavior  of  s^  . 

Plastic  models  in  which  the  yield  condition  is  a 

function  of  ,  Eq  .  (1-4),  can  adequately  represent  this 
*) 

behavior  .  Presently,  the  variable  moduli  model  which  is 
described  by  Eqs.  (111-79)  -  (III-80)  in  unloading  does  not 
adequately  describe  this  behavior  at  very  low  stress  levels. 

f 

It  is  not  clear  that  should  be  the  criterion  upon  which 
the  choice  of  the  proper  G  is  based.  Further  study  may  in- 

.  i 

dicate  that  some  combination  of  and  p  should  be  used  for 
this  purpoi e . 

The  slope  of  the  uniaxial  stress-strain  curve  in  un- 
4 

loading,  j  G^  ,  could  not  be  less  than  unless  G^ 

were  permitted  to  be  negative  for  some  range  of  stresses. 

Thus,  in  order  to  obtain  the  sharp  break  found  in  the  ex¬ 
perimental  curve,  Fig.  (la),  at  low  stress  levels,  a  des¬ 
cription  of  such  that  it  could  be  negative  at  very  low 

s  t  i.*  e  S  S  1  c  V  b  1  S  oppcdTS  tO  ue  promising.  Work  On  tiic  p  F  O  p  e  l* 


It  can  be  shown  that  the  minimum  value  of  the  deviator 
occurs  when  upon  unloading,  the  opposite  face  of  the 
yield  surface  is  reached.  Upon  continued  unloading,  the 
stress  path  is  along  the  yield  surface. 


representation  of  unloading  for  variable  moduli  models  is 
continuing. 

For  the  present  model,  K„.t  must  be  chosen  larger  than 

u 

the  maximum  value  of  K  found  during  loading.  A  second  r  e- 

3 

quirement  is  that  x  must  be  greater  than  the  maximum  value 

•k  ) 

of  .  Since  the  work  is  still  in  the  early  stage  of 

development  and  is  subject  to  revision,  a  detailed  de¬ 
scription  of  the  equations  applicable  to  unloading  will  not 
be  given. 

( e )  Numerical  Results  . 

Typical  results  for  the  variable  moduli  model  of  this 

section  are  shown  in  Figs.  (16)  -  (19).  The  parameters  used 

K1  K2 

in  the  computations  were  V  =  0.30,  —  =  -100,  and  =  4000, 

ON  N 

O  O 

the  same  as  were  used  in  the  previous  strain  model.  In 
addition,  the  values  Yj  =  60  and  y  =  -133.3  were  chosen. 

In  loading,  the  uniaxial  stress-strain  curve,  Fig.  (16),  has 
the  characteristic  reversal  of  curvature  which  is  found  in  ex¬ 
perimental  curves.  Again,  the  point  of  inflection  occurs  at 

a  strain  close  to  4 7-,  a  typical  value.  On  unloading,  with 
KUN 

■  =  30,  the  stress  decreases  shrrplv,  Although  the  un- 
Ko 

loading  portion  appears  to  be  a  straight  line,  the  alone  at 
low  stress  levels  is  in  fact  less  than  half  the  value  at 
high  stress  levels.  Nevertheless,  the  distinct  tail  which 
is  found  on  unloading  experimentally,  Fig.  (la),  does  not 

^  This  corresponds  to  the  requirement  V  >  0. 
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appear,  illustrating  the  inadequacy  of  the  present  unloading 
model  at  these  very  Low  stresses. 

The  plot  of  the  radial  stress  versus  axial  stress  in 
uniaxial  strain  is  shown  in  Fig.  (17).  On  loading,  the 
curve  is  actually  concave  upward  at  very  low  stress  levels. 

At  higher  stresses,  the  curve  is  essentially  a  straight  line. 
The  unloading  radial  stress  is  always  greater  than  the  corre¬ 
sponding  value  in  loading.  The  unloading  curve  is  concave 
downward  and  changes  in  curvature  at  0^  =  0.164,  where  s^ 
changes  sign.  The  plot  has  the  same  general  characteristics 
as  the  experimental  curves,  Fig.  (lb),  except  that  in  the 
experimental  unloading  curve  a ^  drops  off  much  more  sharply 
as  G^  is  brjught  back  to  zero. 

Finally,  the  deviator  s^  is  plotted  versus  the  pressure 
p  in  Fig.  (18)  for  loading  and  unloading  in  uniaxial  strain. 

On  loading,  the  initial  curvature  is  concave  downward.  How¬ 
ever,  the  major  portion  of  the  curve  is  essentially  a  straight 
line.  On  unloading  s^^  is  always  less  than  its  corresponding 
value  in  loading  and  the  curve  is  concave  upward.  At  s^  =  0, 

the  slope  is  continuous,  but  the  curvature  suddenly  increases, 

.  ■ 

reflecting  the  change  in  sign  of  3^  in  Eq .  (III-80)  .  The 
plot  ends  at  »  -p  (so  that  G^  =  0)  with  the  slope  almost 
horizontal . 

The  results  for  the  trlaxial  compression  test  are 
illustrated  for  loading  only  in  Fig.  (19).  The  curves  are 
drawn  for  the  same  parameters  which  were  used  in  the  uniaxial 
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strain  test  and  for  the  same  values  of  the  lateral  stress 

wh'eh  were  used  in  the  strain  variable  moduli  model,  Fig.  (14), 
.  °3  - 

namely  •£—  =  0.04,  0.06,  0.08.  Each  of  the  curves  is  concave 
o 

downward  and  approaches  asymptotically  the  value  (o^  - 
given  by  Eq.  (III-72).  At  a  higher  value  of  the  lateral 
stress,  the  stress  difference  at  failure  increases,  as  does 
the  initial  slope. 

On  the  basis  of  the  present  results,  one  sees  that  the 
theoretical  combined  variable  moduli  material,  when  subjected 
to  two  special  loading  configurations,  namely  the  uniaxial 
strain  and  triaxial  compression  test,  reproduces  all  the 
salient  features  found  experimentally  in  these  tests.  There¬ 
fore,  the  present  model  offers  promise  of  being  able  to  give 
a  reasonable  representation  of  real  soils  in  more  general 
loading  configurations. 

The  tentative  unloading  model  appears  to  mirror  the  uni¬ 
axial  strain  experimental  results  over  most  of  the  stress 
range  of  interest.  Further  work  is  required  to  clarify  the 
material  behavior  at  very  low  stresses  and  to  study  unloading 
in  the  triaxial  compression  test. 

The  ability  of  the  present  model  to  match,  numerically, 
real  soil  data  and  the  process  used  to  determine  the  various 
constants  must  av/ait  the  completion  of  current  investigations, 
both  theoretical  and  experimental. 
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It  is  of  interest  to  explore  in  some  detail  the  relation 
between  the  variable  moduli  models  in  Section  III  and  the 
various  plastic  models  discussed  in  Sections  I  and  II.  This 
discussion  follows  in  Part  (f). 

(f )  Comparison  of  Variable  Moduli  and  Plastic  Models. 

In  the  previous  portion  of  this  section,  sir  rities 
between  the  variable  moduli  models  and  the  plastic  models 
have  been  mentioned.  This  subsection  discusses  several  of 
these  similarities  for  simplified  models  of  these  types.  On 
the  basis  of  the  present  study,  it  appears  that  the  concepts 
of  a  "yield  condition"  and  of  plastic  flow  may  be  contained 
within  the  theory  of  the  variable  moduli  models. 

For  many  materials,  empirical  evidence  suggests  the 
existence  of  states  of  stress  (and/or  strain)  at  which  the 
material  undergoes  continuously  increasing  deformations  with 
little  or  no  increase  in  loading.  This  combination  of 
stresses  at  which  flow  occurs  is  often  called  a  "flow  con¬ 
dition"  or  a  "yield  condition".  When  these  deformations 
become  sufficiently  large  so  that  unacceptable  changes  in 
the  geometry  occur,  this  state  is  called  "failure".  Plastic 
material  models  describe  the  stress  state  at  which  flow 
begins  hy  a  yield  condition  and  the  subsequent  deformations 
by  a  flow  rule.  The  variable  moduli  models  describe  the 
behavior  of  materials  as  they  approach  this  critical  state 
of  stress  (and/or  strain)  as  well  as  their  behavior  at  the 
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state  itself.  For  example,  if  one  were  to  define  an 

"incremental  shear  modulus"  as  the  change  in  shear  stress 

divided  by  the  corresponding  change  in  shear  strain  and 

apply  this  definition  to  the  alluvium-playa  model  in  tri- 

axial  impression,  Fig.  (10),  this  incremental  shear  modulus 

would  be  zero  at  "failure".  Large  increases  in  shear  strain 

would  suddenly  occur  with  no  increase  in  shear  stress  in  this 

situation.  In  the  case  of  either  of  the  variable  moduli 

models,  Figs.  (14),  (19),  the  material  does  not  suddenly  fail, 

but  instead  gets  increasingly  softer  up  to  the  "failure" 

point,  i.e.,  G  decreases  continuously  to  zero.  Examining 

the  behavior  of  the  deviator  s^  for  a  Coulomb  material  in 

uniaxial  strain  in  the  very  low  stress  range  during  unloading 

(this  corresponds  to  the  flat  portion  of  Fig.  2),  one  finds 

that  s^  is  opposite  in  sign  from  0^  and  becomes  smaller  in 

2 

magnitude.  The  strain  deviator,  =  5  ei  f°r  uniaxial  strain, 

dsi 

decreases  monotonica lly .  Thus,  the  ratio  is  negative  in 

this  region.  Although  the  present  unloading  relation,  Eq. 

(III-80),  does  not  adequately  represent  this  behavior  at  low 

stresses,  it  is  felt  that  improved  versions  of  the  mathematical 

*) 

model  in  which  G  on  unloading  may  be  negative  for  certain  com¬ 
binations  of  stresses,  may  provide  a  better  representation. 

An  interesting  illustration  of  the  relation  between 
plastic  and  variable  moduli  models  can  be  obtained  by  a  com¬ 
parison  of  the  simple  Prager-Drucxer  material  and  the  combined 

In  the  present  combined  variable  modulus  model,  the 
shear  modulus  G  is  necessarily  positive  in  both  loading 
and  unloading. 


variable  moduli  material  with  =  K2  =  0.  For  both  models, 
the  Mohr  failure  envelope  is  a  straight  line.  In  the 
Prager-Dr ucker  material,  permissible  status  of  stress  may  be 
defined  in  terms  of  the  yield  condition 


k  +  3ctp 


(See  Eq.  1-1) 


In  the  combined  variable  moduli  material,  permissible 
states  of  stress  are  those  for  which  C  _>  0,  or,  dividing 
Eq.  (III-57 )  by  -y,  >  0 


Equation  (III  —  84)  is  not  an  obvious  requirement  either; 
it  comes  either  from  the  restriction  $  <_  90°  in  plane 
strain,  or  that  the  slope  in  uniaxial  strain  in  un¬ 
loading  be  positive. 


Imm* 


«  * 
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In  uniaxial  strain,  the  requirement  for  initial  soften¬ 
ing,  Eq.  (III-78),  with  =  0  reduces  to 
G  Y, 

2  _£  +  /3  _i  >  ° 

0  Y1 

K 

which  upon  substitution  of  Eq .  (III-83)  and  8  *=  3  r— 

o 

becomes 


a8  < 


/3 


(IT i  —  85) 


(IIX-86) 


The  identical  condition  must  be  satisfied  for  a  Prager- 
Drucker  raaterial  to  yield  in  uniaxial  strain. 

Finally,  examining  the  slope  in  uniaxial  strain  for 
the  combined  mixed  moduli  material,  Eq.  (III-65),  when 
R1  ®  Kj  »  0,  one  finds 


do 


2  Y1 


2  - 


2  Yl, 


77 -  -  [v  G  +  -  ~  K  ]  exp  (—  y.e  )  +  K  (1  ~  “ - ~) 

dGl  3  0  /3  yx  °  /3  1  1  °  /3  Yjl 

The  initial  slope,  Eq .  (III-87)  evaluated  at  =  0,  is 


the  elastic  constrained  modulus  K  +  ~s  G 

o  3  o 


At  large 

strains,  since  y^  <  0,  the  exponential  term  vanishes  and 
the  slope  approaches  asymptotically 


dCj 

de. 


'-YjE  >>  1 


2  Y1 
K  (1  -  —  — ) 

°  /3  Yx 


If  the  and  terms  were  kept  in  the  bulk  modulus 
relation,  the  slope  would  simply  continue  to  increase. 


(111-87) 


(III-K8) 


■> 

{ 
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Tn  the  linear  elastic  Prager-Drucker  material,  the  slope 

in  uniaxial  strain,  [see  Fig.  (2)],  is  the  elastic  slope 
4 

Kq  +  -j  Gq  until  yield,  at  which  point  the  slope  becomes 
discontinuously  the  plastic  slope 


Plastic 


(1  +  2/3  ct)2 
0  (1  +  3a2$) 


The  limiting  value  of  the  variable  moduli  slope 


V1  - 


K  (1  +  2/3  a) 
o 


(111-89) 


(111-90) 


is  not,  however,  the  plastic  slope,  Eq .  (III-89).  The  two 

2 

are  equal  at  the  end  points  a  *-  0  and  a  =  -  of  the  range 

2  e/3 

0  <  a  <  -  ,  but  the  variable  moduli  value  is  slightly  less 

8/3 

tnan  the  plastic  value  elsewhere  in  the  range.  The  value, 

Eq.  (III-90),  is  that  which  would  be  obtained  for  a  hybrid 
material,  i.e.,  one  which  obeyed  a  yield  condition  of  the 
Coulomb  type,  but  obeyed  a  flow  rule  of  the  von  Mises  type,  [2]. 
For  this  material,  even  incremental  dilatancy  effects  (see 
Section  I)  are  suppressed.  It  is  not  surprising  that  in  a 
material  such  as  the  variable  moduli  material  in  which  the 
bulk  modulus  is  a  function  of  the  mean  strain  only  and  in  which 
shear  effects  cannot  cause  an  increase  in  volume,  the  slope 
approaches  that  of  the  hybrid  plastic  material  and  not  that 
of  the  Coulomb  material  itself. 


Finally,  comparing  both  the  triaxial  stress  and  the  uni¬ 


axial  strain  tests  fcr  the  variable  moduli  material  and  the 
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Prager-Drucker  material,  Fig.  (20).  one  sees  that  the  latter 
describes  yield,  while  the  former  describes  material  behavior 
approaching  "yield". 

( 4 )  Closing  Remarks . 

A  theory  of  variable  moduli  materials  has  been  partially 
developed  in  this  section.  The  results  obtained  for  two  of 
the  models  appear  to  essentially  match  those  found  experimentally 
in  both  the  uniaxial  strain  and  triaxial  compression  tests,  at 
least  for  some  soils  in  certain  ranges  of  stress.  Further 
work  must  still  be  done  in  unloading.  Moreover,  as  new  ex¬ 
perimental  results  become  available  for  higher  stress  ranges, 
the  extension  of  the  variable  moduli  models  to  these  higher 
stress  ranges  must  be  examined. 

It  should  be  noted  that  whereas  the  alluvium-playa  model 
discussed  in  Section  II  may  be  considered  to  be  a  relatively 
advanced  plastic  model,  the  variable  moduli  models  presented 
in  this  section  are  a  first  attempt  at  such  a  theory.  Further 
development  will  hopefully  lead  to  variai le  moduli  models  of 
more  geneia.l  applicability. 
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FIG.  5  UNIAXIAL  STR£  iN  -  ELIMINATION  OF  NET  DILATION 


FIG. 6  -  LOAD -UNLOAD -RELOAD  RELATIONS 

IN  MATERIAL  WHICH  SUSTAINS  DILATANCY 
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FIG.  14.  TRIAXIAL  COMPRESSION  TEST 
VARIABLE  MODULI  MATERIAL,  K  =Ko+K,e  +K2e*  ;  G  =  G*  +  G,e  +  G,^ 

(DRAWN  FOR  J/.-0.30,  K|/Ko«-IOO,  K2/K«.-4000,  G,/*.  =  48,  G,/Ko»  -32) 


VARIABLE  MODULI  MIXED  MODEL,  K  =  K„+  K,  e  +  K2e*  6  =G0+X  P  +  X/  'fjl 
(DRAWN  FOR  «0.30,  K,/K,*-IOO,  K„/K<,  =  4000,  X=60,  /  = -133.3, Xuu/K, 
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STRESS  DEVIATOR  VS.  PRESSURE  IN  UNIAXIAL  STRAIN  __ 

VARIABLE  MODULI  MIXED  MODEL,  K  =  K04K,e  +  Kae2jG=G0-+  Y,  p 4- 

(DRAWN  FOR  V0  =  0.30,  KYI*,,  *  -100,  K,/ko  =  4000,  Yi  =60,  >?=-t33.3,  Xoh/ko=30) 
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